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ABSTRACT 

We present theoretical models of black hole (BH) populations in young stellar environments, such as 
starbursts and young star clusters. Using a population synthesis approach we compute the formation 
rates and characteristic properties of single and binary BHs for various representative ages and choices 
of parameters. We find that most of the BHs (typically 80% for an initial 50% binary fraction) are 
single, but with many originating from primordial binaries (which either merged into a single massive 
star or were disrupted following a supernova explosion). A smaller but significant fraction (typically 
20%) of the BHs remain in binary systems. Main-sequence stars are the most frequent BH companions, 
but massive BH-BH binaries are the next most numerous group. The most massive BHs found in 
our simulations reach ~ 80 Mq, and are formed through mergers of massive binary components. If 
formed in a dense star cluster such a massive stellar BH may become the seed for growth to a more 
massive ( "intermediate- mass" ) BH. Although we do not include dynamical interactions, our results 
provide realistic initial conditions for A'^-body simulations of dense star clusters (e.g., globular clusters) 
including primordial BHs. 

Subject headings: binaries: close — black hole physics — gravitational waves — stars: evolution 



1. INTRODUCTION 

Until recently it was believed that BHs existed in two 
separate mass ranges: stellar-mass BHs, with masses 
~ 10 Mq, and supermassive BHs, with masses 10® — 
10^ M0 (see, e.g.. Fryer & Kalogera 2001; Peterson 
2003; Tegmark 2002). However, over the last few 
years, evidence has been mounting for the existence of 
intermediate-mass BHs (Miller & Colbert 2004). Al- 
though the observations remain controversial a variety of 
plausible formation scenarios for intermediate-mass BHs 
(IMBHs) have been proposed (e.g., van der Marel 2003). 
Primordial formation of massive BHs in the early uni- 
verse has been discussed for many years (e.g., Carr 1993; 
Khlopov, Rubin & Sakharov 2002). Later on, when the 
first very massive metal-free (Pop HI) stars form, one 
naturally expects that, given the lack of significant wind 
mass loss, these stars may collapse to form IMBH rem- 
nants with masses up to ~ 10^ — 10^ Mq (e.g., Heger et 
al. 2003). Some of these IMBHs may be in binaries either 
formed through captures in dense environments (Wyithe 
& Loeb 2003) or remaining from primordial Pop HI bi- 
naries (Belczynski, Bulik & Rudak 2004a). Dynamical 
formation processes for IMBHs in Pop II and Pop I star 
clusters have also been the subject of many recent stud- 
ies. The two most promising scenarios involve runaway 
collisions and mergers of massive main-sequence stars 
(Portegies Zwart & McMillan 2002; Giirkan, Freitag, & 
Rasio 2004) and successive mergers of stellar-mass BH 
(Miller & Hamilton 2002). In this paper, we focus on 
the populations of stellar BHs forming out of the most 
massive stars in Pop II or Pop I young stellar environ- 
ments. 

Stellar BHs form through core collapse of massive stars 
with masses > 20 Mq. If N stars form with a standard 
Kroupa (Kroupa & Weidner 2003) initial mass function 



(IMF) between M^in = 0.08 M© and Af^ax = 150 M©, 
we expect A^bh — 5 x 10^'' A^ black holes to have formed 
after ~ 10^ yr, when all stars massive enough to produce 
a BH have evolved. The goal of this paper is to char- 
acterize this initial population of BHs. For a population 
of single stars, this is merely a question of characteriz- 
ing the metallicity-dependent relation between progen- 
itor mass and final BH mass (Fryer & Kalogera 2001; 
Heger et al. 2003). However, most stars, and, especially, 
most massive stars are expected to form in binary sys- 
tems (Duquennoy & Mayor 1991). The BH formation 
processes can be affected significantly by the evolution 
of the progenitor star in a binary, especially if the initial 
orbital separation is < 500 Rq and the star can over- 
flow its Roche lobe as it evolves. In addition, many 
BHs will retain binary companions, which can drasti- 
cally affect their later evolution, their detectability as 
X-ray or gravitational-wave sources, and the way they 
interact with their environment. 

The study we present in this paper is based on a 
population synthesis approach, in which a large num- 
ber of single and binary stars are evolved according to 
parametrized evolutionary prescriptions. Stellar evolu- 
tion is followed starting from a population of zero-age 
main-sequence (ZAMS) stars with specified metallicity. 
Given enough time the stars evolve to form remnants, 
either white dwarfs (not relevant for the massive stars 
considered here), neutron stars, or BHs. Both observa- 
tions and more detailed theoretical calculations are used 
to constrain and parametrize the uncertain stages of evo- 
lution (e.g., supernova explosions and common envelope 
phases). A number of population synthesis codes of vary- 
ing levels of sophistication are currently being used to 
study many astrophysical problems (e.g., Fryer, Woosley, 
& Hartmann 1999; Nelemans, Yungelson & Portegies- 
Zwart 2001; Hurley, Tout & Pols 2002; Pfahl, Rappa- 
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port, & Podsiadlowski 2003; Belczynski, Kalogera & Bu- 
lik 2002, hereinafter BKB02). Here we use the most re- 
cent version of the Star Track population synthesis code 
(BKB02; Belczynski et al. 2004b) to study BH forma- 
tion in young metal-rich (Pop I) and metal-poor (Pop II) 
populations. This work is complementary to the previous 
study of the oldest binary BHs, descendants of metal-free 
Pop III stars, by Belczynski et al. (2004a). 

Our results can be applied directly to the modeling of 
young starburst populations, with ages ~ 10^ — 10^ yr. 
For those systems, we have derived complete synthetic 
sample catalogues of all BHs, single and in binaries, that 
have evolved from the initial burst of star formation. Our 
results can also be applied to the modeling of young star 
clusters (such as the super star clusters so prevalent in 
starbursts; see, e.g., de Grijs et al. 2003), although we 
do not take into account the effects of dynamical interac- 
tions in dense cluster cores (cf. Giirkan, Freitag & Rasio 
2003; Ivanova et al. 2004). 

Characterizing the properties of a primordial BH popu- 
lation is also very important for many theoretical studies 
of old star clusters. Indeed many of these studies attempt 
to bypass the first ^ lO^yr of massive star evolution and 
use initial conditions that already contain BHs. For ex- 
ample, Af-body simulations of globular clusters contain- 
ing primordial BHs have been performed starting with 
a small number of identical single BHs of mass IOMq 
(e.g., Fregeau et al. 2002; Portegies Zwart & McMillan 
2000). One goal of our work is to provide more realistic 
initial conditions for those studies. The further dynam- 
ical evolution of BHs in dense star clusters could play 
a key role in many problems of great current interest, 
such as understanding the formation of IMBHs and ultra- 
luminous X-ray sources (Miller & Hamilton 2002) and 
predicting the merger rate of BH binaries detectable by 
gravitational-wave detectors (Portegies Zwart & McMil- 
lan 2000). 

Our paper is organized as follows. In § 2 we describe 
the evolutionary model used in the population synthesis 
calculations, including details about BH formation. In 
§3 we present the results of our calculations. First we 
present the results for our reference model, then we follow 
with the description of a number of alternative models. 
Finally, in § 4 we summarize our main results. 

2. MODEL DESCRIPTION AND ASSUMPTIONS 

2.1. Population Synthesis Code 

Our calculations are based on a population synthesis 
method. The StarTrack code (BKB02) has recently 
undergone major revisions and updates (Belczynski et 
al. 2004b). These include: a detailed treatment of 
tidal dissipation effects; individual treatments of vari- 
ous Roche lobe overflow (RLOF) phases; full numeri- 
cal orbit evolution with angular momentum losses due 
to magnetic braking, gravitational radiation (GR), mass 
transfer /loss, and tidal interactions, and incorporating 
the stabilizing influence of optically thick winds. 

All stars arc evolved based on the metallicity- 
dependent models of Hurley, Pols, & Tout (2000), with 
the improvements described in BKB02. Each star, ei- 
ther a single or a binary component, is initiated on 
ZAMS, and its evolution is followed from ZAMS through 
a sequence of evolutionary phases: Main Sequence 
(MS), Hertzsprung Gap (HG), Red Giant Branch (RG), 



Core Helium Burning (CHeB), Asymptotic Giant Branch 
(AGB), and for stars stripped off their hydrogen- rich lay- 
ers Helium phases (He). The nuclear evolution of a star 
ends at the formation of a stellar remnant: a white dwarf 
(WD), a neutron star (NS) or a BH. 

We adopt parameters corresponding to the standard 
model of BKB02 (sec § 2 in BKB02) incorporating the 
latest natal kick velocity distribution of Arzoumanian, 
Chernoff & Cordes (2002), and limiting the accretion 
rate onto the NS and the BH at the maximum Eddington 
limit with the rest of the transferred material lost with 
specific orbital angular momentum of the accretor during 
the dynamically stable RLOF events, but allowing for 
hyper-critical accretion during common envelope (CE) 
phases. During dynamically stable phases with non- 
compact accretors, we allow for non-conservative evo- 
lution, with half of the transferred mass lost with the 
specific angular momentum equal to 2'KjA^/P, where P 
is orbital period, A binary separation, and we choose the 
scaling parameter j = 1 for our calculations (Podsiad- 
lowski, Joss & Hsu 1992). In addition to dynamically 
unstable RLOF events we also allow for evolution into 
the CE phase in cases where the trapping radius exceeds 
the Roche lobe radius of the accretor (e.g.. King & Begel- 
man 1999; Ivanova et al. 2003). Phases preceding CE 
are driven by angular momentum losses (as described 
above), while the inspiral during CE is treated in the 
standard manner through a prescribed energy efficiency 
(with Occ X A = 1 in our standard model; for details see 
BKB02). 

2.2. Black Hole Formation 

BHs are formed out of the most massive stars. Using 
the stellar models of Hurley et al. (2000) and Woosley 
(1986) we estimate the time of core collapse of the 
massive star. At that time we know the mass of the 
core (both CO and FeNi core) and the envelope. For 
intermediate-mass stars (^ 20 — 30 Mq for low metal- 
licity models) the FeNi core is collapsed to form a hot 
proto-NS or a low- mass BH, and we use the work of Fryer 
(1999) and Fryer fc Kalogera (2001) to decide how much 
fall back is expected in a given case (based on the mass of 
the CO core). Fall-back material is accreted onto the cen- 
tral objec;t, increasing its mass, while the rest is assumed 
to be ejected in the supernova (SN) explosion. For the 
highest masses (> 30 Mq for low metallicity models), 
the entire star goes into collapse, forming the BH di- 
rectly, with no accompanying SN explosion. The regime 
for direct BH formation may be easily seen from Figure 1 
(pre-coUapse mass equal to the remnant mass). 

The large observed velocities of radio pulsars imply sig- 
nificant asymmetries in SN explosions. Although the un- 
derlying mechanism is not yet understood, it is generally 
accepted that NS can receive substantial kicks at birth 
(~ 100 - 1000 kms-i). Here we adopt the latest NS kick 
velocity distribution of Arzoumanian et al. (2002). Most 
recent observations (e.g., Mirabel & Rodrigues 2003) sug- 
gest that BHs may be formed cither with an accompa- 
nying kick (for smaller mass stellar BHs) or without one 
(for the most massive stellar BHs). This folds naturally 
into our prescription for BH formation. For low-mass, 
fall-back BHs, we expect somewhat attenuated SN ex- 
plosions, and we assume that the kick is smaller as well: 
its magnitude is assumed to be inversely proportional to 
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the mass of fall-back material. For direct BH formation 
(so called silent collapse) there is no explosion, and we 
therefore assume no asymmetry and no natal kick. 

In Figure 1 we show the initial-to-final mass relation for 
single stars for 3 representative metallicities {Z =0.02, 
0.001, 0.0001). Remnants of various types are shown 
with different symbols. We also plot the instantaneous 
mass of the star just prior to the formation of the com- 
pact remnant. 

The top panel of Figure 1 shows the initial-to-final 
mass relation for solar metallicity {Z = 0.02), the case 
described and discussed in detail in BKB02: BHs are 
formed from stars more massive than about 2OM0, with 
a maximum BH mass of IIMq. There is a general rise 
of BH mass with increasing progenitor mass. However, 
since the wind mass-loss rate increases with the mass 
of the star as well (depleting the mass reservoir for BH 
formation) the initial-to-final mass relation flattens out 
for higher progenitor masses. There are two distinctive 
dips, followed by the subsequent flattening of the rela- 
tion. The first (around 25 Mq), less pronounced dip 
corresponds to the point at which single stars are de- 
pleted of their H-rich envelopes, entering the Wolf-Rayet 
stage with enhanced wind mass-loss rates (Hamann & 
Koesterke 1998). Higher mass-loss rates reduce the mass 
of the star and its core, eventually leading to the for- 
mation of lower-mass remnants. The second dip (around 
50 Mq), corresponds to the point where stars reach huni- 
nosities high enough to initiate; the Luminous Blue Vari- 
able (LBV) phase, characterized by extremely high wind 
mass- loss rates (Hurley et al. 2000), and a sudden de- 
crease of BH masses. 

The middle and bottom panels of Figure 1 show 
the initial-to-final mass relations for lower metallicities 
{Z = 0.001 and Z = 0.0001, respectively). In each 
case the relation differs significantly from that obtained 
for higher (solar) metallicity. In particular, there is an 
increase of the maximum BH mass (to ^ 27Mq for 
Z = 0.001 — 0.0001), and the shape of the relation is also 
altered. Most of these changes can be attributed to the 
dependence of the wind mass-loss rates on metallicity. 
Here we have adopted wind mass-loss prescriptions with 
square- root dependence on metallicity (oc ^JZIZq)^ op- 
erating however only during specific evolutionary phases 
(for details see Hurley et al. 2000). The obvious con- 
sequence of smaller mass loss rates for lower metallicity 
is the increase of the stellar mass just prior to the for- 
mation of the remnant, leading directly to higher BH 
masses. In addition, the smaller mass-loss rates imply 
that the stars encounter the LBV dip at a lower mass 
(~ 33M0 for Z = 0.001-0.0001) as they evolve to higher 
luminosities. The opposite is true for the Wolf-Rayet 
dip, as now only the stars with much higher mass lose 
enough material to become naked helium stars (^ 38Mq 
for Z = 0.001 - 0.0001). Note the interesting reposi- 
tioning of the two dips: for solar metallicity the dip at 
high mass corresponds to the LBV phase, while for lower 
metallicity {Z = 0.001 - 0.0001) it corresponds to the 
start of the Wolf-Rayet phase. 

Figure 1 also allows us to see the change in relative 
numbers of BHs formed through fall back or direct pro- 
genitor collapse. This is important for survival of the 
BHs in binary systems, since for fall-back BHs, which 



receive natal kicks, the systems hosting BH progenitors 
may be disrupted by SN explosions. For solar metallic- 
ity many BHs are formed with kicks through fall back, 
which occurs for single stars with initial masses in the 
range 20 - 42 Mq and 48 - 70 Mq. For Z = 0.001, BHs 
receive a kick in the narrower ranges 18 — 25 Mq and 
39 — 54 Mq. For an even lower metallicity oi Z ~ 0.0001, 
only BHs formed from stars in the mass range 18 — 24 Mq 
receive a kick, while all others form silently. 

3. RESULTS 

3.1. Standard Reference Model 

Our reference model starts with 10^ binaries and 10^ 
single stars (initial binary fraction /bi = 50%). The ini- 
tial masses of single stars and binary system primaries 
(more massive) are chosen from a standard initial mass 
function (IMF) with slope —2.35 between 4Mq and the 
maximum mass Mmaxj characteristic of young clusters 
(see Kroupa & Weidner 2003 for detailed discussion). We 
adopt Minax = 150 Mq for the reference model calcula- 
tion. The masses of secondary stars in binary systems 
are sampled from a mass ratio (</, secondary/primary) 
distribution assumed to be constant between and 1 for 
our standard model calculation (but values below the 
hydrogen burning limit at 0.08 Mq are rejected). The 
distribution of initial binary separations is assumed to 
be constant in logarithm between the minimum (such 
that binary components at ZAMS are not in contact) 
and lO^i?0, while for eccentricities we assume a ther- 
mal distribution. In our standard model all stars have 
metallicity Z = 0.001. 

Here we define a BH as a compact object with a mass 
exceeding the maximum NS mass Mmax.NS = 3 Mq (this 
value will be changed later). We examine the BHs at five 
different epochs: 8.7 Myr (corresponding to a turnoff 
mass Mto = 25 Mq), 11.0 Myr (Mto = 20 Mq), 15.8 
Myr (Mto = 15 Mq), 41.7 Myr (Mto = 8 Mq) and 103.8 
Myr (Mto = 5 Mq). The binaries can produce both 
BHs in binaries and single BHs. For the BHs still in 
binaries we list the numbers formed in various binary 
configurations. 

Single BHs originating from binaries are formed 
through several different channels. First, the binary may 
be disrupted in a SN because of mass loss and a possi- 
ble nascent kick. In this case the evolution of the two 
components is followed as two single stars. Single BHs 
formed this way will be listed explicitly under the cat- 
egory "Single: binary disruption". Binary components 
may also merge following close interactions (e.g., CE evo- 
lution, Darwin instability) forming a single object, which 
may later evolve and leave behind a single BH (denoted 
by the category "Single: binary merger"). The merger 
product is assigned a new mass, either the total compo- 
nent masses for mergers involving compact remnants or 
MS stars or their core masses in cases involving giant-like 
stars (with the envelopes assumed lost during the merger 
process). If one component is already a BH, we assume 
that the outcome of the merger is still a BH (formed 
instantaneously, with a new higher mass). In all other 
cases, we assume the formation of a new ZAMS star (i.e., 
"full rejuvenation") and we proceed to evolve the new 
single star. If the merger product is massive enough a 
(single) BH may form. Since all merging binary compo- 
nents are already evolved, the evolution leading to the 
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formation of a BH in this case may in reality be some- 
what shorter, and thus our predicted numbers for this 
category are only lower limits. Most of the merger prod- 
ucts massive enough to produce BHs are found to have 
come from pairs of massive HG-MS, He-MS, and MS-MS 
systems. 

In Table 1 wc list the numbers of BHs found both as 
single objects and as members of binary systems. Bi- 
nary BHs are most frequently found either in BH-MS 
or BH-BH systems. Shortly after the initial starburst, 
MS companions naturally dominate (they are the most 
probable BH companions at early times). However, as 
the age increases these massive MS companions them- 
selves evolve, collapse to BHs and some eventually form 
BH-BH systems (dominating at later times). The sec- 
ond longest lived phase in massive star evolution (after 
the MS) corresponds to core helium burning. We note 
a similar trend for BH-CHeB as for BH-MS binaries. 
Although naturally less numerous overall, the number 
of BH-CHeB binaries peaks as soon as these systems 
appear but then tends to decrease with time. As the 
less massive progenitors end their lives and start forming 
remnants, we observe increasing numbers of BH-NS and 
later BH-WD systems. 

There is an overall decrease in the number of binary 
BHs with time. This can be understood as follows. Most 
massive stars have finished their evolution and formed 
BHs at early times (< 10 Myr). Many of the BH binaries 
with unevolved MS companions either merge in RLOF 
events or are disrupted when the secondary star under- 
goes a SN explosion. At later times, some of the tightest 
binaries may also merge due to orbital decay driven by 
GR emission. The opposite trend — increase with time — 
is seen in the total number of single BHs. After the first 
^ 10 Myr most of the single BH progenitors have fin- 
ished their evolution and have formed BHs; therefore the 
later increase is connected to binary evolution. The bi- 
nary mergers and disrupted components continue adding 
to the single BH population at later times. For an initial 
50% binary fraction we predict that single BHs eventu- 
ally dominate by about an order of magnitude over BHs 
found in binaries. 

There is an interesting direct relationship between the 
relative numbers of binary and single BHs for any spe- 
cific starburst and its initial (primordial) binary frac- 
tion. Etc. Let us define the ratio of binary to single 
BHs (as a function of the initial binary fraction /bi) as 
-R(/bi) = (100% iVbi„)/7Vsin, where A^biii and iY^in denote 
the total number of BHs found in binaries and as single 
objects, respectively. Obviously, this ratio is also a func- 
tion of age, but, for definiteness, let us pick a specific age 
of t = 11 Myr. For an initial population consisting en- 
tirely of binaries (/bi = 100% from Table 1 we can read 
that the total number of binary BHs is 27729 while the 
number of single BH is 64271 (all from binary disrup- 
tions and mergers), so the ratio i?(/bi = 100%) = 0.37. 
This is obviously an upper limit. Any contribution from 
single stars in the initial population will lower this value. 

3.2. Parameter Study and Normalization 

To assess the effects of various model assumptions on 
BH formation, we repeat our calculations for a number 
of different models, each differing from our standard ref- 
erence model in the value of one parameter or one as- 



sumption. The specific models wc chose arc described in 
Table 2. 

Normalization: In addition to a brief description of 
each model. Table 2 also lists for each simulation the 
total initial mass of single and binary stars. All stars 
are assumed to form in an instantaneous burst of star 
formation. Stars can form from the hydrogen burning 
limit up to the maximum mass Mmax characterizing a 
given system. We adopt the three-component, power- 
law IMF of Kroupa, Tout, & Gilmore (1993) with slope 
ai = —1.3 within the initial mass range 0.08 — 0.5 Mq, 
a2 = —2.2 for stars within 0.5 — 1.0 Mq, and = —2.35 
within 1.0 Mq — Mmax- This IMF is easily integrated to 
find the total mass contained in single and binary stars 
for any adopted ai,a2 values. The particular choice of 
low-mass end slope of the IMF (Q!i,a2) usually (with 
the exception of Model B) does not change our results. 
Low-mass stars do not contribute to the BH populations. 
However, as most of the initial stellar mass is contained 
in low-mass stars, a small change in the IMF slope at 
the low-mass end can significantly change the normaliza- 
tion. Similarly, our results can easily be generalized to 
other primordial binary fractions (/bi) by simply weigh- 
ing differently the results obtained for single stars and 
for binaries. 

Summary of results: In Table 3 we summarize the state 
of each simulation at 11 Myr (Afto = 20 M©) for the 
models presented in Table 2. In the following, we de- 
scribe in more detail how the changes made in each model 
affect the various BH populations. 

In Model B we create the binary by independently (1) 
drawing the primary from 4 — 150Mq according to an 
IMF exponent -2.35 and (2) drawing the secondary from 
0.08 — I5OM0 via the Kroupa broken power-law IMF. 
This initialization process results in the ratio q of com- 
ponent masses peaked at low q values. For this model, al- 
most all BHs in binaries are found with MS companions. 
In addition, relative to the reference model, fewer single 
BHs are produced from disrupted binaries. Since the ini- 
tial mass ratio is on average smaller in Model B, massive 
stars, the progenitors of BHs, are formed preferentially 
with low-mass MS companions. If a BH progenitor does 
not destroy the binary, then the low mass companion re- 
mains for a long time on MS (large number of BH MS 
systems) and it does not initiate either RLOF nor ex- 
plode as SN (small numbers of single BHs). In contrast 
to the standard model, there are almost no BH BH sys- 
tems, as they originate typically from binary progenitors 
with relatively high mass ratios that are almost absent 
in Model B. 

For Model CI we use small metalUcity Z = 0.0001 

(characteristic of old Pop II stars metallicity) while for 
Model C2 we choose high (solar) metallicity of Z = 0.02. 
We find that the total number of BHs formed can vary 
strongly with the initial chemical abundance. While sin- 
gle BHs are not greatly affected, the population of bi- 
nary BHs changes drastically with metallicity. With in- 
creasing metallicity the evolution of the star is altered; 
most notably, the stellar wind mass loss rate increases, 
causing the star's mass to drop more rapidly. Eventu- 
ally, some stars either form lower mass BHs or instead 
of forming BHs they end up as NSs. Also higher wind 
mass loss rates (Model C2) and their associated angu- 
lar momentum losses leave binaries wider, and therefore 
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i) easier to disrupt (see the increase of single BHs from 
binary disruption), ii) making component interactions 
harder (smaller chance of RLOF; see the smaller merger 
BH number). To summarize, then, the BH population 
formed from binary stars depends sensitively on the char- 
acteristic metallicity of the primordial population. 

In Model D we decrease the CE efficiency to 0.1, which 
means that only 10% of binary orbital energy may be 
used for envelope ejection. This leads to either merger 
of binary components in CE phase, or results in a much 
tighter orbit of post-CE binary, if it survives. There is 
no significant change in the number of BHs in Model D 
as compared to the reference model. Only the close bina- 
ries with rather extreme mass ratios undergo a dynam- 
ically unstable RLOF, and thus their evolution may be 
altered in Model D. This obviously depletes number of 
BH MS stars (as they tend to have rather small mass 
ratios, otherwise a MS star would have already evolved 
to a BH), while leaves the number of BH-BH systems 
unchanged (formed at early times out of comparable 
mass binary components). Further orbital shrinkage dur- 
ing the dynamically unstable RLOF events obviously in- 
creases the number of binary mergers (especially in the 
second RLOF episode), but the merger products do not 
have enough time to form BHs after 11 Myrs. ^ Fi- 
nally the binaries which evolved through one or two CE 
events are tighter and are harder to disrupt in SN ex- 
plosions (explaining the slight decrease in the number of 
BHs from binary disruptions). 

In Model E, we apply for all BHs, independent of their 
mass and formation scenario, full nascent kicks following 
the distribution of Arzoumanian et al. (2002). By con- 
trast, in our standard model, we either have used smaller 
kicks or no kicks at all for BH formation. Therefore, 
we expect (and observe) a drastic decrease in number of 
BHs formed in binary systems, since full kicks disrupt 
all wide and also many close binaries. Further, since sys- 
tems that would otherwise be binary BHs are disrupted, 
we produce many more single BHs than in our reference 
model. Although the nascent kicks play an important 
role in the modeling of BH binary populations, we note 
that the most massive BHs almost certainly form silently 
(Mirabel & Rodrigues 2003), without nascent kicks, so 
the results of Model E should be treated as an extreme 
case. 

In Model F, for high-mass stars (potential BH pro- 
genitors), we use a steeper IMF exponent {a^ = —2.7), 
similar to the one observed in field populations. As ex- 
pected, all the numbers for BHs found in binaries and as 
single objects are decreased. Although we do not expect 
the initial starburst IMF to be as steep as the one ob- 
served in the field (Kroupa & Weidner 2003) , it is worth 
noting the factor of 2 decrease in BH numbers for this 
choice. 

Since smaller, less massive clusters form stars only up 
to a smaller maximum mass, in Models Gl and 02, we 
decrease Mmax to 50 and 100 Mq, respectively. A de- 
crease in the maximum stellar mass in a given simulation 
shifts more stars to lower mass, depleting the number of 
potential BH progenitors. Therefore there is a general 
decrease of BH numbers found in Models G. Since more 

^ Note that at later epoch the number of binary mergers has 
increased by factor of ~ 3 (see Table 4). 



BH progenitors are removed from initial population in 
Model Gl, it shows the largest change as compared to 
the standard model calculation. However, due to a steep 
IMF for massive BH progenitors, the high mass end miss- 
ing in Models Gl and G2 does not contain many stars 
and the change is not large. 

In Model H we show results for a lowered maximum 
mass for NS formation, which may be as low as 1.8 — 
2.3 Mq (see, e.g., Akmal, Pandharipande, & Ravenhall 
1998). In this model we adopt Mmax.NS = 2M0, and 
all compact objects over that mass are assumed to be 
BHs. As expected, this model has slightly more single 
and binary BHs. Clearly the maximum NS mass is not 
significant compared to other model uncertainties. 

Finally, under "Model I" (which is not really a different 
model), we give the numbers of BHs with mass greater 
than lOAfQ found in our reference model. Since most of 
the single stars and non-interacting binary components 
form BHs with M > 10 Mq (see Fig. 1), there is only a 
slight decrease in BH numbers. 

In Table 4 we present the subpopulations of single and 
binary BHs at 103.8 Myr (Mto = 5 M©). The trends are 
very similar to these observed at the earlier epoch (Table 
3), however, with the increased contribution of remnant 
BH binaries (BH-NS and BH-WD). 

3.3. Orbital Periods of Black Hole Binaries 

Period distribution for the standard model: In Figure 2 
we plot the distribution of orbital periods of BH binaries 
for our standard model calculation at f = 11 Myrs after 
starburst. This figure also shows the contributions from 
the BH MS and BH BH subpopulations. 

Orbital periods of binaries hosting BHs are found in 
a very wide range Porb ~ 0.1 — 10^ days. The distri- 
bution is characterized by two distinctive peaks, smaller 
at Porb ~ 10 days, and larger at Porb ~ 10^ days, with 
the underpopulated region (Porb ^ 10'^ days) in between. 
The clear distinction between short and long period bina- 
ries is caused by progenitor system RLOF history or lack 
of thereof. There is a bifurcation period at Pbur ~ 10^ 
days: below which systems go at least through one RLOF 
interaction leading most frequently to the orbital shrink- 
age, while the wider systems never interact, but only lose 
material in stellar winds, and their orbits expand. The 
specific value of Pbur is different for each binary, and de- 
pends on the maximum radii the components may reach 
in their evolution (set by initial mass and metallicity) . If 
at any moment of the binary evolution radius of either 
component increases over its Roche lobe radius (set by 
the period and mass ratio) then RLOF starts. Therefore, 
the Pbur is a function of a number of parameters describ- 
ing the binary and its components and it depends also 
on the orbit evolution (e.g., tidal interactions), and it is 
found within a range: Pbur ~ 100 — 10^ days for binary 
BH progenitors. Since the tighter binaries evolve through 
one or more RLOF interactions, many of them merge, 
and single objects are formed, causing the depletion of 
the short-period binary BH population (smaller peak). 
Although wide binaries avoid mergers, they are much 
more prone to the disruption in SN explosions. However, 
since most of BHs are formed silently or with rather small 
kicks, the depletion is not large, and the peak at long pe- 
riods is substantially higher than for tighter binaries. 

The BH-BH binaries evolve only through GR. Their 
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orbits decay slowly, and for the tightest binaries the two 
components can merge and produce a single BH (within 
a Hubble time). 

Systems containing a non-remnant companion (e.g., 
BH-MS, BH-CHeB), given enough time, will populate 
other BH subpopulations. In particular, if the compan- 
ion is massive enough to form a second BH in the system, 
and the binary survives potential interactions and SN 
explosion, a BH-BH binary will be formed. It is worth 
noting that the majority of the binaries are quite "soft," 
with long periods Porb ^ 100 days, so that, if placed 
in a dense stellar environment (e.g., a globular cluster) 
they will be easily disrupted, further enhancing the sin- 
gle BH population. On the other hand, a smaller but 
significant number of binaries are hard and, through dy- 
namical interactions with other stars, may even become 
harder (evolve to shorter periods). The overall shape of 
the distribution remains mostly unchanged at different 
times (listed in Table 1). However, as discussed in §3.2, 
the contribution of the BH-BH binaries increases, while 
the BH-MS systems are relatively less frequent at the 
later times. 

Dependence of period distribution on model assump- 
tions: In Figure 3 we present period distributions of BH 
binaries for all our alternative evolutionary models. 

In most cases (e.g. Models B, CI, D, F, Gl, G2 
and H), the general two-peaked shape of our reference 
model distribution is preserved. For some of these mod- 
els position and/or relative hight of the two peaks may 
change slightly, for reasons that can be easily explained 
in terms of the different model assumptions. For exam- 
ple in Model D, we find that the number of tight binaries 
in comparison with long-period systems is much smaller 
than for the reference model. Moreover, the gap between 
the two peaks is more pronounced and wider than for the 
reference model. Model D simulations used smaller CE 
efficiency, affecting tight interacting binaries by shifting 
them to smaller orbital periods (wider gap) or leading to 
component merger (smaller peak at low periods). 

For a few models, the changes in the distribution are 
much more pronounced. In Models C2 (high metallicity) 
and I (high mass BHs only) the distribution is peaked at 
large orbital periods, while the short period BH binaries, 
are relatively very infrequent. For example, the Model I 
evolution (which favors only the highest mass progeni- 
tors) leads to a selection and survival of only the widest 
binaries. If the high mass stars, progenitors of the most 
massive BHs, are placed on the tight orbits, they merge 
in the first RLOF interaction, and the formation of the 
tight BH system is aborted. In contrast with previous 
two, in Model E, we find that almost all of the BH bina- 
ries have short orbital periods. With full kicks applied 
to all BHs, almost all wide systems are disrupted. 

3.4. Black Hole Masses 

Mass distribution for the standard model: In Figure 4 
wc plot the mass distribution of the entire BH popula- 
tion formed in oiu' standard Model A. Additionally, we 
plot separately the single and binary BH subpopulations. 
The distribution shows three well defined peaks, first at 
Mbh ~ 6 - 8 M0, second at Mbh ~ 10 - 16 Mg, and 
third at Mbh ~ 22 - 26 M©, and then it steeply falls off 
with the increasing BH mass. The distribution may be 
divided into two major contributions, one from the popu- 



lation of single BHs, and the other from BHs in binaries. 
We sec that single BHs dominate the population (see also 
§3.2) and basically set the shape of the overall mass dis- 
tribution. A majority of the single BHs originate either 
from single star progenitors or from disruptions of wide 
(non-interacting) binary stars in SN explosions. There- 
fore, the shape of distribution for single BHs is mainly 
determined by a combination of the initial-to-final mass 
relation for single BHs (Fig. 1) with the IMF. And, in 
fact, in Figure 1 (middle panel) we can see the pile up 
of BH remnants corresponding to the three peaks in the 
BH mass distribution of Figure 4. Stars over ^ 50 Mq 
form 10 — 16 M© BHs corresponding to the second, the 
largest, peak in Figure 4. Stars with initial masses of 
25 — 35 Mq form ~ 25 Mq BHs, providing a majority 
of the third peak in the BH mass distribution. Finally, 
the stars with initial masses of 40 — 50 Mq tend to form 
BHs with masses of 7 Mq, the first peak of the Fig- 
ure 4 distribution. The initial-to-final mass relation lead- 
ing directly to the shape of the BH mass distribution is 
described in detail in §3.1. The binary evolution may 
increase or decrease the mass reservoir for BH formation 
through RLOF interactions between system components. 
However, with one or two exceptions, the overall shape of 
the mass distribution for single and binary BHs is rather 
similar. Specifically, the most apparent change is that 
for binary BHs the first narrow peak for BH masses of 
6 — 8 Mq is missing. This peak for single BHs comes 
from the turnover in the initial-to-final mass relation, re- 
sulting in a pile up of BHs in this specific mass range. 
This characteristic feature of the initial-to-final mass re- 
lation corresponds to a very sharp transition in single 
star evolution, from H-rich to naked helium stars, which 
is caused by wind mass loss and the effective envelope re- 
moval for single stars above a certain initial mass. Since, 
for binary stars, removal of the envelope is allowed not 
only through stellar winds but also through RLOF in- 
teractions, it is allowed for the entire mass range, and 
therefore there is no sharp transition between evolution 
of H-rich and helium reach stars, and the aforementioned 
peak disappears. 

One more important effect of binary evolution on BH 
final mass is illustrated in Figure 5, where we show sin- 
gle BHs up to the highest formed BH mass. Single BHs 
form either from single progenitors, from components of 
disrupted binaries, or through binary mergers. Clearly 
we can see that single star progenitors form BHs only up 
to ^ 30 Mq (as expected from Fig. 1). Slightly higher 
BH masses 40 Mq) are obtained through binary dis- 
ruption, as some of the progenitors may have been reju- 
venated in RLOF event to a higher mass in the preced- 
ing binary cn'olution. All of the most massive BHs, up 
to ^ 80 Mq, are formed through binary mergers. Most 
mergers are formed out of HG MS pairs, due to the very 
rapid expansion in the Hertzsprung Gap, leading very 
often to dynamically unstable RLOF (in which case we 
always assume a merger; see Ivanova & Taam 2004). In 
case of such a merger, the final mass of the newly formed 
single star is the sum of the MS star mass and the core 
mass of the HG star (the envelope of the HG star, con- 
taining > 50% of the entire star mass, is assumed to be 
lost in the merger process). Since we allow for significant 
mass loss in the merger process, the final merger remnant 
mass may be underestimated. Had we assumed that all 
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mass in the system remains after merger, some of BHs 
produced through this channel would have even larger 
masses, > 100 Mq. These objects have high enough 
masses to be classified observationally as "intermediate- 
mass" BHs. However, it is important to note that they 
are formed through ordinary binary star evolution, with- 
out the need for additional dynamical interaction pro- 
cesses. 

Effects of model assumptions on mass distribution: In 
Figures 6 and 7 we present the BH mass distributions 
for all the models in our parameter study. The overall 
distribution for the entire BH population is not greatly 
affected by different choices of parameter values, with the 
exception of metallicity (see Fig. 6). This is easily under- 
stood, as the highest mass BHs are formed only at the low 
metallicity (Models CI and A) and the lightest BHs are 
formed at the high metallicity (Model C2) as discussed 
in § 3.1. Single BH masses (Fig. 7) in many models reach 
very high values around 80 Mq, and this is a robust 
result of our calculations for a number of different evo- 
lutionary and initial conditions. The highest maximum 
BH masses are found for lowest metallicity environments, 
in a larger systems (with high Mmax), and for binaries 
formed with flat mass ratio distributions, quite indepen- 
dent of other evolutionary parameters. Only in a few 
models, with high metallicity ( C2), uncorrelated initial 
binary component masses (B), and low Mmax (Gl, G2), 
does the maximum BH mass stay below ^ 50 Mq. 

GR emission: BH-BH binaries (rather than double 
NSs), are probably the best candidates for detection 
by ground-based interferometers (Lipunov, Postnov, & 
Prokhorov 1997; BuHk & Belczynski 2003). BH BH sys- 
tems therefore are an important candidate for present 
projects to detect astrophysical GR sources (LIGO, 
VIRGO). The properties of BH-BH binaries at differ- 
ent metallicities formed within much larger stellar sys- 
tems with continuous star formation (e.g., Galactic disk) 
were extensively studied previously (Bulik & Belczynski 
2003; Bulik, Belczynski & Rudak 2004a; Bulik, Gondek- 
Rosinska & Belczynski 2004b). We find that the proper- 
ties of BH-BH binaries in starbursts are not too different 
from those found in the previous studies. Most of BH- 
BH systems are characterized by rather equal masses, 
with a mass ratio distribution peaking at q ^ 0.8 — 1.0 
(for comparison see low-metallicity models of Bulik et al. 
2004b). For most models only a small fraction (~ few per 
cent; e.g., 2% for Model A) of the BH-BH systems are 
tight enough to merge within a Hubble time and produce 
observable GR signals. However, for models which tend 
to produce tighter BH-BH systems (D, E) the fraction 
can be significantly higher (~ 10-20 %). 

4. SUMMARY 

We have calculated the evolution of the massive stars 
found in young stellar environments, such as starbursts. 
The final products of massive star evolution, single and 
binary BHs, were then studied. All our results were dis- 
cussed taking into account the many model uncertainties. 
A number of alternative calculations with varied initial 
conditions and evolutionary parameters were performed 
and presented. We also supplied the necessary data to 



normalize our results to any given total mass of a star- 
burst galaxy or star cluster, with arbitrary choice of ini- 
tial binary fraction. The calibrated results may then be 
used as part of initial conditions for realistic A''-body sim- 
ulations of dense stellar clusters that include primordial 
BHs. 

Soon after the initial star formation burst, most BHs 
are found as single objects, although a significant frac- 
tion of BHs are also found in binaries. A number of single 
BHs are formed as the end product of binary evolution, 
either through binary disruption following a SN explo- 
sion or through a merger following a dynamically unsta- 
ble RLOF episode. The most common binary BHs are 
BH-MS and BH-BH systems. The period distribution of 
binaries containing BHs is usually bimodal, with a ma- 
jority of systems in the long period peak {P ~ 10* — 10^ 
days). These wide binaries would be "soft" if placed in 
a dense stellar environment (e.g., a globular cluster) and 
they would then be disrupted following any strong in- 
teraction with another passing star or binary, thereby 
further enhancing the population of single BHs. The re- 
maining, short-period BH binaries {P ^ 1 — 100 days) 
would instead undergo hardening and evolve, over many 
relaxation times, to produce a population of very com- 
pact binaries that could eventually merge through GR 
emission. 

The typical BH masses are found to be within the range 
7 — 25 Mq, both for single and binary BHs. However, 
the single BHs formed through binary mergers can reach 
masses as high as ~ 80 M©. Since most mergers are 
assumed in our models to be accompanied by significant 
mass loss, BHs formed through binary evolution (without 
any dynamical interactions) could in principle reach even 
higher masses, up to ~ 100 Mq (in the absence of signif- 
icant merger-induced mass loss). This result has many 
important implications. First, some ultra-luminous X- 
ray sources might be explained by a ^ 100 Mq stellar 
BH accreting from a lower-mass companion. This would 
require the capture of a new companion (most likely 
through an exchange interaction with another binary), 
but no dynamics would be involved in the BH formation 
(cf. Kalogera, King, & Rasio 2004). Second, these most 
massive stellar BHs may act as seeds for the formation 
of true IMBHs (with masses > 1000 Mq) that could re- 
side at the centers of some dense star clusters (Gebhardt, 
Rich, & Ho 2002; Gerssen et al. 2002; Miller & Hamil- 
ton 2002). Third, a broader mass range for the tightest 
BH-BH binaries (possibly undergoing further hardening 
through dynamical interactions in a dense star c;lustca') 
will modify predictions for the gravitational- wave signals 
detectable by laser-interferometer instruments such as 
LIGO and VIRGO (Flanagan & Hughes 1998). 
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Type^ 


8.7 Myr 
Mto = 25 Mq 


11.0 Myr 
Mto = 20 Mq 


15.8 Myr 
Mto = 15 M0 


41.7 Myr 

Mto = 8 Mq 


103.8 Myi 

Mto = 5 M, 


BH-MS 


17315 


16207 


12215 


6571 


4004 


BH-HG 


22 


16 


14 


16 


7 


BH-RG 











1 





BH-CHeB 


1254 


1029 


675 


262 


155 


BH-AGB 


16 


13 


27 


9 


16 


BH-He 


167 


102 


60 








BH-WD 











1 


1075 


BH-NS 


69 


364 


760 


913 


880 


BH-BH 


9261 


9998 


10022 


10010 


9996 


Total in binaries: 


28104 


27729 


23773 


17783 


16133 


Single: binary disruption 


24093 


43909 


55649 


60262 


61021 


Single: binary merger 


10128 


20362 


33825 


65236 


66148 


Single progenitor 


77580 


108080 


120100 


120100 


120100 


Total single: 


111801 


172351 


209574 


245598 


247269 



''Black holes in binary systems are listed aecordiiig to their companion types: MS — main sequence, HG — Hertzsprung Gap, RG — rcg 
giant, CHeB — core He burning, AGB — asymptotic giant branch. He — helium star, WD — white dwarf, NS — neutron star, BH — black hole. 
Single black holes formed from components of disrupted binaries are listed under "Single: binary disruption." Single black holes formed 
from binary merger products are under "Single: binary merger." Single black holes that are remnants of single stars are listed under "Single 
progenitor." 



Table 2. Population Synthesis Model Assumptions 



Mass [Mq] in Mass [Mq] in 



Model 


Description* 


Single Stars 


Binaries 


A 


standard model described in §3.1 


3.8 X 10'^ 


5.9 X lO'^ 


B 


uncorrclated binary component masses 


3.8 X 10^ 


7.6 X 10^ 


Cl-2 


metaUicity Z = 0.0001,0.02 


3.8 X 10^ 


5.9 X 10'^ 


D 


QfCE X A = 0.1 


3.8 X 10^ 


5.9 X 10'' 


E 


full kicks for BHs 


3.8 X 10^ 


5.9 X 10' 


F 


steeper IMF: as = -2.7 


6.0 X 10^ 


9.5 X 10' 


Gl 


lower maximum mass: Mmax = 50Mq 


3.7 X 10'^ 


5.8 X 10' 


G2 


lower maximum mass: Afmax = IOOMq 


3.8 X 10'^ 


5.9 X 10' 


H 


Mmax.NS = 2M0 


3.8 X 10'^ 


5.9 X 10' 


fb 


BHs more massive than 10 Mq 


3.8 X lO'^ 


5.9 X 10' 



''Details of model assumptions are given in § 3.2 and § 3.3. 

''Model I is shown only to give the numbers of BHs (formed in the standard Model A) with mass greater than IOMq. 
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Table 3. Very Young Black Hole Populations - Parameter Study'' 



Type 


A 


B 


Gl 


G2 


D 


E 


F 


Gl 


G2 


H 


I 


Binaries: 
























BH-MS 


16207 


32464 


25527 


2831 


12683 


3176 


9656 


13434 


15798 


17566 


11024 


BH-HG 


16 


I 


26 


7 


25 


g 


10 


17 


28 


16 


10 


BH-RG 



































RH-GHpR 


1029 


16 


1571 


200 


967 


58 


616 


776 


1032 


1047 


894 


J_>J.i. 


13 


Q 


20 




25 


Q 


10 


13 


20 


13 


12 


BH-He 


102 


3 


204 


71 


26 


76 


64 


102 


116 


106 


22 


BH-WD 



































BH-NS 


364 


2 


545 


94 


99 


106 


180 


305 


364 


55 





BH-BH 


9998 


26 


20914 


2743 


10129 


196 


4784 


3127 


7496 


10326 


7674 


Total: 


27729 


32512 


48807 


5953 


23954 


3618 


15320 


17774 


24854 


29129 


19636 


Single: 
























binary disruption 


43909 


18929 


21037 


38110 


34806 


79260 


25696 


35858 


42022 


47714 


10863 


binary merger 


20362 


7402 


16790 


33809 


20084 


19830 


10552 


10837 


18692 


20793 


18030 


Single progenitor 


108080 


108080 


109755 


122470 


108080 


108080 


63090 


83695 


102255 


108080 


87630 


Total: 


172351 


134411 


147582 


194389 


162970 


207170 


99338 


130390 


162969 


176587 


116523 



'AH numbers correspond to an age of 11 Myrs {Mto = 2OM0) 
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Table 4. Young Black Hole Populations - Parameter Study^ 



Type 


A 


B 


CI 


C2 


D 


E 


F 


Gl 


G2 


H 


I 


BiricLriGSi 
























BH-MS 


4004 


31762 


5952 


888 


3052 


766 


2520 


3394 


3998 


4506 


2836 


BH-HG 


7 


2 


14 





11 





6 


5 


14 


7 


5 


BH-RG 











1 








2 














BH-CHgB 


155 


50 


238 


55 


153 


2 


96 


104 


164 


155 


151 


BH-AGB 


16 


3 


22 


6 


18 





8 


16 


14 


16 


15 


BH-He 





























1 





BH-WD 


1075 


205 


1470 


459 


1064 


10 


612 


779 


966 


1087 


1033 


BH-NS 


880 


13 


1705 


160 


292 


356 


512 


763 


980 


610 


56 


BH-BH 


9996 


26 


20915 


2743 


9745 


202 


4790 


3130 


7494 


10434 


7674 


Total: 


16133 


32061 


30316 


4312 


14335 


1336 


8546 


8191 


1630 


16816 


11770 


Single: 
























binary disruption 


61021 


23660 


44889 


40146 


51089 


88394 


36806 


51708 


58766 


88195 


15167 


binary merger 


66148 


19012 


54904 


88961 


74952 


64172 


43794 


57573 


65004 


68834 


49515 


Single progenitor 


120100 


120100 


147065 


127835 


120100 


120100 


93080 


124755 


142765 


148465 


87630 


Total: 


247269 


162772 


246858 


256942 


246141 


272666 


173680 


234036 


266535 


305494 


152312 



*A11 numbers correspond to an age of 103.8 Myrs [Mto = 5Mq) 
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Fig. 1. — Initial-to-final mass relation for diflferent metallicities. Remnants of different type: white dwarfs (WD), neutron stars (NS) 
and black holes (BH) are marked. 
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Fig. 2. — Period distribution of BH binaries (solid line) for our standard model after 11 Myr. Two major contributing system types 
arc shown separately: BH-MS binaries (dashed line) and BH— BH binaries (dotted line). Normalized to total number of BHs (single and 
binaries), binwidth: tenth of decade. 
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Fig. 3. — Period distribution of BH binaries (solid line) for different models, all at 11 Myr (dashed line; BH-MS binaries; dotted line: 
BH-BH binaries). 
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Fig. 4. — Mass distribution of BHs at 11 Myr (solid line) for standard model ( A). Single BHs are shown with dotted line and BHs in 
binaries with dashed line. Normalized to total number of BHs (single and binaries); binwidth: 0.5Mq. 
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Fig. 5. — Mass distribution of various l<inds of single BHs at 11 Myr in standard model ( A). The dotted line shows BHs originating 
from primordial single stars; the dashed line represents single BHs from disrupted binaries; the solid line is for single BHs that are the 
remnants of merged binaries. Normalized to total number of BHs (single and binaries); binwidth: 2.0Mq. 
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Fig. 6. — Mass distribution of BHs at 11 Myr for different models. Conventions are as in Fig. 4. Note that for Models CI and C2 the 
vertical scale differs from the rest of the panels. 
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Mass distribution of various kinds of single BHs at 11 Myr for different models. Conventions are as in Fig. 5. 



